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Introduction SCA in Practice Estimating Mixing Matrix 
Radio frequency interference (RFI) is a problem for microwave remote sensing of Earth. eas S,(t) 


Scatter plot of 
coefficients c, 


Although frequency allocations are set aside for passive sensing, RFI can still degrade X,(t) 
measurement quality. In some cases radiometer bandwidth exceeds allocated Pi 
spectrum to reduce measurement uncertainty or spectrum allocations are shared, 
forcing microwave radiometer to co-exist with terrestrial sources. Low level RFI is 
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particularly detrimental as it can be concealed as natural variability leading to flawed X,(t) ™ C,,(t) C,,(t) 
scientific results. RFI detection algorithms have been developed to address the * Each source is transformed to a sparse Algorithms 
problem. Research into other algorithms is needed to improve upon the sensitivity of representation in another domain. * Dictionary * The mixing matrix A is estimated through the scatter plot of the coefficients c,. 


existing detection algorithms to various types of RFI. The Sparse Component Analysis 
(SCA) has been investigated to determine its sensitivity to continuous wave (CW) RFI. 


Sparse Component Analysis (SCA) 


SCA is a blind source separation method which seeks to extract N unknown sources 


* Dictionary learning 
* Structured dictionary 
« Joint sparse representation 
* Global optimization 
* Greedy algorithm 
* Mixing matrix estimation 
* Global clustering 
* Local scatter plots 
* Separation 
+ Binary masking 


Mixing matrix A is estimated using the sparse 
coefficients obtained in the transformation. 


The sources are separated in the transformed 
domain using the coefficients of observations and 
the estimated mixing matrix A. 


The sources are finally reconstructed in the time 


A global clustering algorithm, weighted histogram, is used to estimate the 
directions of the columns of the mixing matrix. 


Separation and Reconstruction 


from P observations where P< N. The sources need to have disjoint supports. domain. * Binary masking is used to separate the sources. 
x(t) = As(t), t=1,..,T n . : . ve. 5! : 
ee os © 7 © RFI Detection * Amask is derived by setting the indices of the sources which correspond to the 
observations mixing matrix sources Xy = Qy4Ny + ayyny highest amplitude in observations given the mixing matrix with a value of one. 
(known) (unknown) — (unknown) Ho: Xy = d24Ny + ag2ny * Since the signals are disjoint, only the active sources are assigned a coefficient 
No RFI * The detection criterion is the median of the absolute value of the reconstructed value, while the other sources are masked with zeroes. 
Ny sources, Z in time. 
Ge) = M12 ean ny * The output of SCA are three reconstructed sources in time, 5; (t), $(t), 3(t) where * Inthe last step, the acquired source signals can be reconstructed back to their 
xy M21 [22 423/\ | _ f=4, 2: 3... Nv original domain by using x = ®c,. 
Hy: i rae one isi * The median of the absolute value of each reconstructed source 
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Scatter plot of the mixtures x, 
and Xp 


Use least squares to estimate columns 
of mixing matrix, A 


<x(t),An> ~ 
pe ars An 
Least squares: §,(t) = ||n| 
0 otherwise 


In practice the source signals are not 
disjoint in time. Columns of A not easily 
determined from scatter plot 
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and INRy = INR, +3dB 


Sparse Domain Transform 


Represent each source by a linear combination of a few elementary signals (atoms): 
S(t) = Dhar s(x (6) 


Measure of sparsity 


1/t 
lleslle = [ks Ies(k)I*) t > 0 quantifies the sparsity of c, 


num of nonzero coefficients t = 0 
Ideal t for sparsity is 0 


Dictionary 


Dictionary: set of elementary signals known as atoms ® 

Consider dictionaries that span C’, K > T: redundant dictionary, infinite 
representation 

Once the dictionary is determined the signal representation can be determined in 
the transformed domain using either a global optimization technique or greedy 
algorithms. 

Let K, denote the number of atoms in the dictionary,®. If the signal in time domain 
is x and the signal coefficients in the transformed domain is c,, then 


coofficients of observations 
0 


x= @c, 2° 


where x = [x(1),...,x(t)]” and c, = [c,(1), ....¢, CK)". ‘1 ae 


cy = arg min |lcx|lq : 
. x=@c, * of sources mixture 
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Conclusion 


* Monte Carlo simulations with 1000 time samples are used, along with the 


structured dictionary and Orthogonal matching pursuit (OMP) algorithm for joint 
sparse representation, the weighted histogram for matrix estimation and binary 
masking for source separation. 


+ Figure shows the performance results of SCA for detection of CW RFI with NR, 


ranging from -15 dB to 2.5 dB. 


* The results show perfect to near perfect detection for INRs greater than -12.5 dB 


and very good detection at -12.5 dB and -15 dB. 


+ Results show that detection works for relatively large INRs for CW RFI. 
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